Bioelectromagnetic calculations in anisotropic volume conductors 
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1 Introduction 

Bioelectromagnetic measurements are widely used 
to obtain structural and functional information non- 
invasively on the human brain (EEG, MEG) and the 
human heart (ECG, MCG). Often it is important to 
locate the origin of the measured signal. The localiza¬ 
tion of the source is called the inverse problem. Solv¬ 
ing the inverse problem is often done iteratively. This 
solving method requires accurate calculation of the 
electromagnetic fields produced by the chosen source 
model. This is called the forward problem of bioelec¬ 
tromagnetism. 

Human tissue is often modeled as being homogen¬ 
eous and isotropic in its electrical properties. How¬ 
ever, this is not always a realistic assumption. In some 
tissue types, the conductivity is strongly direction- 
dependent. Examples of this are the myocardium, 
the brain, and ordinary muscle tissue. In the myocar¬ 
dium, as well as in ordinary muscles, the anisotropy is 
due to the cell structure of the muscles, which favors 
the electric conduction into the direction of the cells 
and hinders the conduction in the perpendicular direc¬ 
tions. The human brain, on the other hand, is mainly 
constructed from nerve cells, whose axons are good 
electrical conductors. The axons tend to form clearly 
oriented bundles, as is seen in Fig. 1. The bundles 
introduce complex anisotropy into the brain [1]. 

The purpose of this study is twofold: 1) To imple¬ 
ment efficient and accurate numerical methods for the 
anisotropic forward problem, and 2) to analyze the 
effects that anisotropy of the volume conductor in¬ 
troduces into the electromagnetic fields of a dipole 
source, the emphasis being on the second point. 

2 Methods 

2.1 The finite element method 

A finite element method (FEM) was implemented in 
order to solve the electric potential produced by a 
current dipole. The major advantage of FEM over 
the standard boundary element method (BEM) is that 


FEM can take into account anisotropic material prop¬ 
erties. In FEM, the volume conductor is divided into 
three-dimensional elements and for every element, an 
individual conductivity tensor can be defined. The 
FEM in this work was implemented using the Galer- 
kin method [2], in which the basis functions are also 
used as test functions in the integral equations. 

The electric potential u produced by a current source 
J in a conducting body f l satisfies the equation 

V • <tV« = V ■ J, (1) 

where o is the conductivity (tensor). A boundary con¬ 
dition for the potential is found by restraining all the 
volume current inside the volume H 

(aVu(r)) ■ n= 0, redfl, (2) 

where n is the outward unit normal vector of the sur¬ 
face dil. The solution of Eq. 1 satisfies also the ex¬ 
pression 

f ^(r)V • (<rVu) tKl = f ^(r)(V • J) dLl, (3) 
Jn Jn 

where Eq. 1 is simply multiplied by the test function 
and integrated over the volume O. 

Next, orthogonal basis functions are introduced. 
The basis functions are defined using the volume ele¬ 
ment mesh. To the zth vertex of the mesh, a single 
basis function is attached. The function has 
the value one at the zth vertex and the value zero at 
all other vertices. The basis functions have a limited 
support; a single basis function has values other than 
zero only in the volume of the elements that are adja¬ 
cent to the particular vertex. In this work, piecewise 
linear basis functions were used. 

The basis functions can be used to approximate any 
well-behaving function in the volume LI. For the elec¬ 
tric potential u the approximation yields 

N 

u = Y^u{n)ipi, (4) 

where N is the number of vertices, u is the approx¬ 
imation of u, and r, is the position vector of the zth 
vertex. 






Figure 1: Lateral aspect of the right hemisphere 
of the brain with arrows pointing out a nerve fiber 
bundle [3]. 


Figure 2: Isocontour lines of the electric potential on 
the surface of the sphere: A) the analytic, and B) the 
numerical solution. The values are millivolts and the 
depth of the source is 20 mm. 

from the Biot-Savart equation 


An interesting point in finite element modeling of the 
forward problem is the form of the source term. We 
used the current dipole as a source model, which can 
easily be included into the FEM equations. The cur¬ 
rent dipole at point ft is defined by J = Q8(r— ft), 
where Q is the dipole moment. Now the source term 
of Eq. 3 can be written 


[ ipi(r)V ■ J dLl = -[ Vtpi(r)-Jdn 

J n Jn 

= - [ Vipi(r) ■ QS(r— ft) dfl 

J n 

= ~Q ■ V<pi(ft), (5) 


where the first equality follows from vector manipula¬ 
tion and the assumption that there is no source current 
on the surface of the conductor. The third equality fol¬ 
lows from the properties of the Dirac delta function. 
The above derivation allows us finally to write the 
weak form equation 


- / Vtpi ■ (crV^Pj) dn=-Q- Vpj(ft), 

r i Jn 

( 6 ) 

which must hold for every basis function <$, Vi = 
1,..., i\T, The equation can also be expressed in mat¬ 
rix form 

Au = b, (7) 

where the matrix A and the vectors u and b are ob¬ 
tained straight from the previous equation. The finite 
element solver implemented in this work constructs 
Eq. 7 and solves it for the potential u. 


2.2 Calculation of the magnetic field 

The magnetic field outside the volume ft, originating 
from a current source in the volume, can be calculated 


B(r) = ^ [ Jtot(r') x -^3 dv', (8) 
47T J v ' r— r'\° 

where J to t is the total volume current, and r and r' 
are the field point vector and the integration variable, 
respectively. The magnetic field produced by a source 
current J can be solved when the total current dens¬ 
ity is divided in two terms: J to t = J — erVu. Since 
the potential u is known as the finite element solution, 
Eq. 8 can be integrated for the magnetic field. In this 
work, the numerical integration was performed using 
a Gaussian quadrature rule for every element separ¬ 
ately, and summing the single element contributions 
together. 

2.3 Volume conductor modeling 

Two different volume conductor shapes were used in 
the study, namely the sphere and the hexahedron. The 
sphere was used to represent the brain, which was 
assumed homogeneous except an anisotropic arc ap¬ 
proximating a nerve fiber bundle. The radius of the 
arc was 25 mm and its cross-sectional area ca. 1 cm 2 . 
The sphere was divided into 31930 linear tetrahed¬ 
ral elements with 6063 nodal points. The element 
mesh density was highest near the north pole of the 
sphere where the anisotropy was included, and the 
mesh density decreased toward the south pole of the 
sphere. The included anisotropic arc consisted of 161 
elements; the anisotropy ratio of 1:10 was used. 

The hexahedral volume conductor was used to model 
the thorax. The volume was also meshed with alto¬ 
gether 31944 linear tetrahedral elements using 6348 
nodal points. The top layer (thickness 8 mm) of the 
hexahedral volume was defined to be tangentially an¬ 
isotropic in the ratio of 1:14. This construction was 




used to approximate the skeletal muscle surrounding 
the human thorax. 

The current source in this work was the point-like cur¬ 
rent dipole. In both volume conductors, a few differ¬ 
ent source locations were studied. In the sphere, the 
dipoles were tangential and at the depth of 13-33 mm, 
the radius of the sphere being 80 mm. In the hexa¬ 
hedron, the dipoles were also tangential, and at the 
distance of 10-36 mm from the top surface, the depth 
of the volume being 240 mm (other dimensions were 
400x400 mm). 

3 Results 

3.1 Validation of the numerical methods 

The accuracy of the field calculations was validated 
by comparing the numerical results to the analytic 
solutions available in the homogeneous conductors. 
The results in the hexahedral volume were compared 
with the solutions in the semi-infinite space. 

The comparisons yielded good correspondence of the 
field patterns measured with the correlation coeffi¬ 
cient p, and fairly good correspondence in the amp¬ 
litudes of the fields. The results of the comparisons 
are collected into the Table 1. An example of the cal¬ 
culated fields is shown in the Fig. 2. 

Table 1: Comparison between numerical and analytic 
results. Shown are the correlation coefficient p and 
the change of the maximal amplitudes of the fields 
in percent in the numerical solution compared to the 
analytic solution. 



Sphere 

Hexahedron 


P 

Ampl. / % 

p Ampl. / % 

Electric 

0.98 

90-190 

0.96 80-250 

Magnetic 

0.99 

90-100 

1.00 100-110 


As Table 1 shows, the numerical methods were quite 
accurate. The values of the correlation coefficient 
show good agreement in the solutions, but in the amp¬ 
litudes of the fields more variability is found. Prob¬ 
ably, using a denser element mesh would improve the 
results also in this regard. However, we believe that 
the amplitude comparisons between numerical results 
in the isotropic and in the partly anisotropic media are 
still meaningful, since the solutions are found with the 
same method for both media. 


3.2 Anisotropic arc in the sphere 

The effect of anisotropy on the bioelectromagnetic 
fields was studied by comparing numerically calcu¬ 
lated results in isotropic and in partly anisotropic me¬ 
dia. The results show that the anisotropy has principal 
effect on the fields originating from a source near or 
inside the anisotropic region. Also, the changes due 
to the anisotropy seem to be larger in the electric po¬ 
tential than in the magnetic field. These effects can be 
seen in the Table 2, where data in the sphere model are 
collected. 

Table 2: Changes in the electric potential (E) and 
in the magnetic field (M) due to added anisotropy 
in the sphere. Correlation coefficient and amplitude 
changes are divided into different categories depend¬ 
ing on the distance between the source of the fields 
and the anisotropic region. The label ’near’ means 
that the source is situated in an anisotropic element 
or in an element next to an anisotropic element, and 
’further’ labels all other sources. 


Correlation coefficient | Amplitude / % 



near 

further 

near 

further 

E 

0.78-0.98 

1.00 

100-290 

100-110 

M 

0.94-1.00 

1.00 

90-170 

90-100 


As Table 2 shows, the anisotropy has almost no ef¬ 
fects on the fields when the current source is loc¬ 
ated outside the anisotropic region. However, when 
the source is close to the anisotropy, some deviations 
compared to totally isotropic case are present. Also, 
the amplitudes of the fields in the isotropic and an¬ 
isotropic case are very close to each other when the 
source is not close to the anisotropy. In the other case, 
rather large differences can be found in the fields. 

3.3 Anisotropic layer in the hexahedron 

The results show that an anisotropic layer on top of 
a hexahedral volume conductor changes the electro¬ 
magnetic fields originating from a source situated be¬ 
low the layer. The amplitude of the fields on the sur¬ 
face of the volume decreases from the influence of the 
anisotropic layer. For electric potential, the maximum 
amplitudes diminish to 30-70% of the amplitudes of 
the isotropic case. For magnetic field, the changes are 
less dramatic, the amplitudes decreasing only to 80- 
90%. 

Also the spatial distribution of the electric potential 
changes substantially, the correlation coefficient be- 





Figure 3: The electric potential on the surface of the 
hexahedron: A) isotropic, andB) anisotropic medium. 

ing in the range of 0.85-0.97. The anisotropic layer 
stretches the potential in the direction of the best con¬ 
ductivity (see Figure 3). This effect explains the rel¬ 
atively low values of the correlation coefficient. The 
magnetic field pattern remains practically unchanged 
by the anisotropy as shown by the value of the correl¬ 
ation coefficient, which in this case is 1.00. 

The behavior of the electric potential in this case is 
interesting regarding the localization of the source 
based on electric measurements. The stretching of the 
field pattern increases the distance between the min¬ 
imum and maximum of the field. The field pattern 
could be then falsely assumed to come from a source 
lying deeper than it actually is, because in an isotropic 
medium the distance between the extrema of the field 
is proportional to the depth of the source. Also, the es¬ 
timated magnitude of the source would be incorrect. 

4 Discussion 

In this article, numerical methods for calculating the 
bioelectromagnetic forward problem in an anisotropic 
volume conductor are presented. To our knowledge, 
there are very few groups that have published biomag- 
netic calculations in an anisotropic medium [4], [5]. 
The volume conductor modeling in this study is very 
simplistic but nevertheless the models and the type of 
anisotropy used are closely related to the real anatom¬ 
ical structures present in the human body. 

The results of the study show that the electric poten¬ 
tial is generally more affected by anisotropic struc¬ 
tures than the corresponding magnetic field. In the 
source localization problem, this property of the mag¬ 
netic field may prove to be an advantage over elec¬ 
tric measurements at least in cardiography where both 
methods are able to detect the signals of the same 
sources. 

It was also found that an anisotropic subvolume in an 
otherwise isotropic medium affects the electromag¬ 
netic fields if the source of the fields is inside the 


anisotropic region, or in contact with it. When the 
source is not in connection with the anisotropic re¬ 
gion, even if it is between the electrodes and the an¬ 
isotropy, there is almost no change in the fields. 
Finally, the results suggest that an anisotropic layer 
on top of a volume conductor affects quite signific¬ 
antly the electric potential making the source local¬ 
ization possibly erroneous if anisotropy is not taken 
into account. The magnetic measurements are likely 
to perform better also in this special case. 

We conclude that the anisotropy of the brain should 
be taken into account if one tries to solve the inverse 
problem as accurate as possible. Also, in the case 
of cardiography, the anisotropy of the skeletal muscle 
layer should be considered. It must be noted that the 
present study is carried out with a very simple set of 
models of the real biological geometry. In our opin¬ 
ion, however, the conclusions of this paper give the 
basic idea of the effects of the volume conductor an¬ 
isotropy. 
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